import nibabel as nib
import matplotlib.pyplot as plt
import numpy as np

#加载.nii.gz文件
img=nib.load('1.nii.gz')
image=img.get_fdata()#获取图像数据
header=img.header#获取图像头信息
num_slices=image.shape[2]#获取切片数量
print(num_slices)

# for i in range(0,num_slices,20):
#     slices = image[:,:,i]
#     plt.imshow(slices,cmap='gray')
#     plt.title(f'Slice {i}') # 设置切片标题
#     plt.show()  # 显示所有子图

fig, axes = plt.subplots(nrows=1, ncols=5, figsize=(60, 60))  # 创建子图
for i in range(5):
    slices = image[:, :, i]  # 获取第i个切片
    axes[i].imshow(slices, cmap='gray')  # 在第i个子图中显示切片
    axes[i].set_title(f'Slice {i}')  # 设置切片标题

plt.tight_layout()  # 调整子图布局以防止重叠
plt.show()  # 显示所有子图
